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Abstract 

We show how to include a chemical potential \x in perfect lattice ac- 
tions. It turns out that the standard procedure of multiplying the quark 
fields ^ at Euclidean time t by exp(±//i), respectively, is perfect. As 
an example, the case of free fermions with chemical potential is worked 
out explicitly. Even after truncation, cut-off effects in the pressure and 
the baryon density are small. Using a (quasi-)perfect action, numerical 
QCD simulations for non-zero chemical potential become more powerful, 
because coarse lattices are sufficient for extracting continuum physics. 
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Understanding strongly interacting matter at finite baryon density is a long- 
standing and challenging problem, motivated for instance by relativistic heavy ion 
collision and by the physics of neutron stars. The standard procedure to formulate 
lattice QCD at a finite chemical potential fi includes a factor exp(±/x) in the time-like 
link variables |J. As a consequence, the Euclidean action is complex, the Boltzmann 
factor cannot be interpreted as a probability, and standard Monte Carlo techniques 
fail. 

The usual method to handle a chemical potential is to simulate at fi = 0, and 
to include the baryon number term by some re-weighting technique in measured 
observables |J. However, this method is tractable only on small physical volumes 
V, for a recent review see Ref. ||. The essential numerical problem is to measure 
exponentially suppressed observables, like the partition function ratio Z(/x)/Z(0) ~ 
exp(— (5V[f{pL) — /(0)]). Here (3 is the inverse temperature and /(/i) is the free 
energy density. In numerical simulations the above ratio arises as an average over 
many positive and negative contributions. Hence its accurate determination requires 
tremendous statistics. An improved lattice action can not directly solve this sign 
problem, but it would help because it suppresses the artifacts due to the finite lattice 
spacing. 

As a particularly troublesome effect caused by lattice artifacts, there is an upper 
limit for the possible fermion number density on the lattice. The value of this limit 
depends on the lattice action. It can be understood from the maximal coupling 
distance in the time-like direction, which restricts the exponential growth of the 
leading term at \i — > oo. An improved action can weaken this unphysical saturation 
effect, and push it to a larger chemical potential. 

As a further problem, chiral symmetry is already restored — and a non-zero 
quark number sets in — as /i reaches half the pion mass . Generally the quenched 
approximation has been blamed for that threshold ||, but — at least with respect 
to the quark number — it persists in simulations using the Glasgow method ||. 
It could be a further manifestation of lattice artifacts, which might cause a spu- 
rious scattering of the eigenvalues @. This problem further motivates the use of 
improved lattice actions at finite baryon density, since they render the observables 
more continuum like than a standard lattice formulation. 

Perfect lattice actions provide a possibility to completely eliminate cut-off effects. 
In practice, approximately perfect actions have been constructed, and they have led 
to a significant reduction of finite lattice spacing artifacts ||, |9], pT0|, fTlf] . When applied 
to a perfect action at fi = 0, the standard procedure of multiplying the fermion fields 
^ at time t by exp(±/it), respectively, turns out to be quantum perfect in the fully 
interacting theory. We use the example of a free quark to demonstrate that — even 
after truncation — the resulting lattice action approximates continuum physics very 
well. [] On coarse lattices, the scaling behavior of observables is strongly improved 

1 Recently another paper has appeared, where a chemical potential was included in a (classically) 
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compared to the standard action and other proposals. For discussions of the free 
standard lattice fermions at /i > 0, see the second Ref. in |l| and Refs. 0, |12 |. 



Let us consider a block factor n renormalization group transformation, which 
maps a fine lattice theory of link variables U and fermion fields \I/ onto a coarse 
lattice theory with fields U', \t>'. The actions S[U, f] and S'[U', f ] on the 
fine and coarse lattice at \x = are related by 

exp(-S"[E/', tf*]) = J VUV^V^exp[-S[U,^,^]-T[U',^'^',U^^}) , (1) 

where T characterizes the renormalization group transformation. It must be chosen 
such that the partition function remains invariant, Z' = Z. Starting on a very fine 
lattice with quark mass m/(nN) and inverse temperature f3nN, and performing N 
renormalization group steps with block factor n, leads in the limit nN —> oo to a 
perfect action S*[U, ^, for quarks of mass m at inverse temperature f3. 

Now we address the question how to incorporate the chemical potential. In the 
continuum this can be done simply by a substitution in momentum space, 

(Jfe) -»■ ^(k, h - in) , * (Jfe) -» fc 4 - i^)- (2) 

In perturbation theory, it is equivalent to leaving the fermionic fields unchanged and 
replacing 

— > A; 4 + i/x (3) 

in the free propagator instead. A priori it is not obvious how to transfer this substi- 
tution to a standard lattice formulation. An early guess replaced sin fc 4 — > sin k^ + ifi 
in the fermionic propagator, but this formulation does not have a proper continuum 
limit. Instead it turned out that one should use exactly the same substitution (0) 
resp. (^|) also on the lattice to incorporate \i in the actions resp. propagators of 
Wilson or staggered fermions 

Hence, on the initial fine lattice a chemical potential fi can be introduced by 
replacing V(x,t) and \I>(f,t) (t e [0,/?]) in the Wilson action S[U, by 

t) = exp(//t)tf (£, t), ^(x, f) = exp(-/it)*(£, t). (4) 

The chemical potential appears as a purely imaginary constant Abelian gauge poten- 
tial A* = z/i. Thus, the fields ^(x, t) and ^(x, t) can be viewed as being parallel 
transported to t = 0. This ensures covariance under the corresponding Abelian 
gauge transformations. In the presence of a chemical potential, the renormalization 
group transformation must be modified if one wants to manifestly preserve this ad- 
ditional symmetry. This is achieved naturally by using T[U', fin ty' ^ n ^} 
with 

^'(f , f) = exp{fmt')^'{x', t'), " n tt'(f, t') = exp{-fint')^'{x', t'). (5) 
perfect action of the 2d 0(3) model p3| . 
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Note that on the coarse lattice t' G [0,/3/n]. Performing the corresponding renormali- 
zation group transformation results in the action S'[U', fJ,n ty'^ n ^']. Indeed, this is 
exactly what one obtains by applying the standard procedure to include a chemical 
potential [in directly in the coarse lattice action S'. Starting with a small chemical 
potential fj,/(nN), and again iterating the blocking procedure N times, yields the 
perfect action S'*[[/, A ' \]/] with chemical potential /i in the limit nN — > oo. Hence, 
performing the standard procedure of including /x in a perfect action — constructed 
at /i = — yields a perfect action at arbitrary fi. We emphasize that this argument 
applies to the fully interacting quantum theory. 



As a special case, in perturbation theory one may send the blocking factor n — > 
oo. Then we don't need to iterate any more; the perfect action is obtained for 

0, 



because one does 



N = 1. We call this technique "blocking from the continuum" 
not start the blocking process from a fine lattice but directly from the continuum 
theory. From that procedure it is particularly evident that the continuum relation 
(0) is inherited without alteration by the perfect lattice propagator. This agrees 
with the result that we derived above using a finite blocking factor. 

For a free fermion with mass m at p, — 0, a perfect action has been constructed 

d 4 k 



in Ref . M 



S*[V,V] = 



b (2tt) 4 

n 2 (A; + 27r/) 



1 

- + -, 

m a 



(6) 



Here, B =] — ir, 7r] 4 is the Brillouin zone and a is a parameter in the renormalization 
group transformation. At a = oo and m = 0, the perfect action is chirally invariant, 
and hence — in agreement with the Nielsen-Ninomiya theorem — the action is 



nonlocal JT^|. At finite a, on the other hand, chiral symmetry is explicitly broken in 
the action — though still present in the observables |j| — and the action becomes 
local. In particular, for a = (e m — m — l)/m 2 , locality - - in the sense of an 
exponential decay of the couplings — is optimal. 

Applying the standard method of including the chemical potential, the perfect 
action at \i > takes the form ^]. In momentum space we obtain 



d A k = 



b (2vr) z 



m(-k)^- 1 (k,k i + i^{k). 



(7) 



Let us consider the cut-off effects in the pressure p and the baryon number density 
riB on the lattice. In the continuum, for massless quarks at zero temperature, these 
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quantities are given by 

/z 4 2/i 3 
P = "7. 9 , n B 



Qtt 2 1 9tt 2 

On the lattice they take the form 

r d 4 k 

p = t — rr \ogdetA(k, fc 4 ) — logdetACk, k 4 + i/j) 

JB (27TJ 4 L 

1 d 

The lattice propagator for a perfect action is given in eq. (|), but one may also 
insert other fermionic lattice propagators known from the literature, such as the 
Wilson fermion, staggered fermions or the so-called D234 action |Tj|. The latter is 



an extension of the Wilson fermion in the spirit of Symanzik's improvement program: 
the artifacts due to the finite lattice spacing are cancelled in the leading order by 
additional couplings along the axes. 

In order to make the perfect action applicable, its couplings have to be truncated 
to a short range. Hence a fair comparison to other lattice fermions can only deal 
with a truncated perfect fermion. In Ref. |T| we performed such a truncation for 



the perfect fermion action (||) — with the parameter a optimizing locality — to the 
couplings in a unit hypercube. Figure 1 illustrates the cut-off effects in the ratio 
p/fj, 4 as a function of \l in lattice units for the zoo of lattice fermions mentioned 
above. In the continuum limit /i — > 0, they all approach the correct value l/6n 2 . At 
finite /i, the Wilson action suffers from a severe cut-off dependence. The behavior 
is somewhat better for staggered fermions, but the artifacts are still quite bad. The 
D234 action is successful at small fi, but around \i = 1.5 it collapses completely. The 
action of the "hypercube fermion" (truncated perfect fermion), however, remains 
close to the continuum value of the considered ratio over a wide range of chemical 
potentials. For example, the considered ratio deviates by less than 50% from the 
continuum value up to /i = 0.8 for the Wilson and the staggered fermion, up to 
fi = 1.6 for the D234 action, and up to fi = 8.2 for the hypercube fermion. We note 
that an alternative truncated fixed point fermion, which has been constructed by 



the Bern/Boulder group |p7| , shows a similar scaling quality; it stays even a little 



closer to the continuum value than the "hypercube fermion" in Figure 1. 

Figure 2 illustrates the behavior of the baryon number density % divided by 
/i 3 . Again we compare the Wilson fermion, the staggered fermion, the D234 action 
and the truncated perfect action. We see that the behavior of our second scaling 
quantity, ns//i 3 , is qualitatively very similar to p/fi 4 . We summarize the outcome 
of our comparison again: the Wilson and staggered fermion are plagued by large 
artifacts even at small fi, the D234 action is good on fine lattices but disastrous on 
very coarse ones, and the "hypercube fermion" stays by far closest to the continuum 
value. As a further example, the baryon number susceptibility xb = dns/dfi leads 
to a scaling quantity Xb/^ 2 , which again behaves similarly for the lattice fermions 
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Figure 1: The ratio p/fi 4 as a function of the chemical potential ji in lattice units. 
The solid line is the continuum value I/Gti 2 . The dashed lines corresponds to the 
Wilson action (r = 1), the staggered fermion action, the D234 action and the trun- 
cated perfect action ( "hypercube fermion"). 



considered above. This is also in qualitative agreement with comparisons done for 
other quantities before: in Ref. [|TTJ we compared another thermodynamic scaling 
ratio, p/T 4 at // = 0, as well as the dispersion relations, and we found in both cases 
qualitatively the same behavior as in Figures 1 and 2. 

Figure 3 illustrates the saturation effect at large chemical potential for various 
lattice actions. At large /i, most lattice actions lead to an asymptotic fermion 
number density rif — 32. That number corresponds to an occupation of each lattice 
point by 16 flavors with spin up and down. This also includes the doubler fermions, 
which are not suppressed any longer at sufficiently large fi. The Wilson action with 
Wilson parameter r = 1 is special, because in that case the time-like doublers are 
completely removed from the spectrum (their mass diverges) [Tj|. As a consequence, 
the asymptotic value for nf is reduced to 16. For the staggered fermions the doub- 
ling is reduced by construction, hence their saturation level is only 8. It should 
be noted that a non-truncated perfect action does not suffer from any saturation 
artifacts; it behaves like the continuum action. This is due to the presence of 
couplings over infinite distances in Euclidean time, even though these couplings are 
exponentially small. In general, the maximal saturation level is proportional to the 
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Figure 2: The ratio ng/fi 3 as a function of the chemical potential // in lattice units. 
The solid line is the continuum result 2/9tc 2 . The dashed lines corresponds to the 
Wilson action (r = 1), the staggered fermion action, the D234 action and the trun- 
cated perfect action ( "hypercube fermion"). 



largest time separation with non-zero couplings. f\ Since the hypercube fermion has 
been truncated at time separation 1, it saturates at n/ = 32. From this point of 
view, the D234 action could, in principle, reach n/ = 64 because it couples over time 
separation 2. However, similar to the effect described above for the Wilson action 
with r = 1, the actual saturation level is only half of that. Q 

Figure 3 further reveals that, as /x grows, the saturation occurs soon for the 
r = 1 Wilson fermion, and at about the same point for the D234 fermion. Staggered 
fermions saturate at even smaller /i. The saturation can be delayed, however, for 
Wilson fermions with r close to (but not equal to) 1, which generates a very heavy 
class of doublers. Then nj reaches a first plateau at 16, performs another jump 
as catches up with the heaviest doubler mass, and finally saturates at nf = 32. 
Also for the hypercube fermion the saturation is significantly delayed, and again we 
recognize a jump as fi exceeds the heaviest doubler mass, so that the final saturation 

2 This provides an argument that any (exactly) perfect action (in more than one dimension) 
must couple over infinite distances. 

3 What happens technically is that the apparently dominant term at ^ — > oo cancels for the 
Wilson fermion at r — 1 and for the D234 action. 
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Figure 3: The fermion number density nj as a function of /i for various lattice 
fermions. The Wilson action with r ^ 1, the D234 action and the hypercube fermion 
saturate atrif = 32, whereas the Wilson fermion with r = 1 (the staggered fermion) 
only reaches nj = 16 (rif = 8). Compared to the usual fermion actions, the sa- 
turation is delayed significantly for the hypercube fermion. For the perfect fermion 
action, n/ is not bounded. 

We have shown how to incorporate a chemical potential in a perfect lattice action 
for full QCD. In an application to free fermions we observed a strongly improved 
scaling behavior at finite /i, even after truncation. We further confirm the scenario 
that the unphysical saturation effect is tamed for a (truncated) perfect action. It is 
straightforward to extend this study to more general cases, involving for instance a 
finite quark mass Q, a finite temperature or an anisotropic lattice. The improvement 
permits lattice simulations for QCD at finite baryon density to be performed on much 
coarser lattices than it is the case for a standard lattice action. Although there is 
also the major sign problem due to the complex action, such an improvement is 
important. 



4 A large fermion mass could be of interest along the lines of the static quark method []19] 
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